Orientational symmetry-breaking correlations in square lattice t-J model 
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We study a tendency to orientational symmetry breaking of the square lattice t-J model by 
applying the exact diagonalization technique. We introduce a small external anisotropy to t and 
J, and calculate the induced anisotropy of hole-hole correlations as a response. It is found that 
the response is strongly enhanced for particular band parameters, meaning that the system has a 
strong tendency to orientational symmetry breaking. The analysis of the momentum distribution 
function indicates that the correlation develops when the Fermi surface is close to (tt, 0) and (0, tt). 
Some properties of high-T c cuprates previously attributed to stripes can also be understood from 
orientational symmetry breaking alone. 
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High-Tc cuprate superconductors are realized by a car- 
rier doping into Mott insulators, where strong onsite 
Coulomb repulsion prohibits double occupation and one 
electron resides at every sites. Such strong electron- 
electron correlations are widely recognized to be crucial 
to the physics of high-T c cuprates. 

Minimal models for high-T c cuprates are believed to 
be the two-dimensional (2D) t-J and Hubbard models on 
a square lattice. Cuprate superconductors have d-wave 
gap symmetry, [l| which was successfully reproduced in 
these models. The models, however, exhibit several or- 
dering tendencies competing with d-wave superconduc- 
tivity: antiferromagnetism,|2, 01 spin-charge stripes, Q| 
staggered flux,[{| and d-density wave.0 These possible 
competing orders were recognized as playing crucial roles 
in high-T c cuprates and discussed in relation to pseu- 
dogap phenomena, vortex structures, and charge 
ordering. [IfJ 

Recently, through microscopic analysises of the t-J [Tlj 
and Hubbard [l2l flii. IbH Hfij models another competing 
order was proposed, the d-wave Fermi surface deforma- 
tion (dFSD): the Fermi surface (FS) expands along the 
k x direction and shrinks along the k y direction or vice 
versa. This order is generated by forward scattering of 
electrons close to the FS near (tt, 0) and (0, tt), and hence 
should be distinguished from all the competing orders 
mentioned above, which require large momentum trans- 
fer with q k (tt, tt). 

Some authors termed the spontaneous dFSD as Pomer- 
anchuk instability, referring to a stability criterion for 
normal Fermi liquids by Pomeranchuk. [16| Note, how- 
ever, that the dFSD turned out to occur usually through 
a first order transition at low temperature, namely with- 
out breaking Pomeranchuk's stability condition. [171 lisj 

The dFSD competes with other instabilities. In partic- 
ular, it can be overwhelmed by d-wave singlet pairing'. [ill 
ITflj However, it was pointed out that even if the sponta- 
neous dFSD was prohibited, the correlations survived po) 
and made a system sensitive to an external anisotropy, 
which led to a noticeable dFSD. [Til l2lj This scenario 



was invoked for high-T c cuprates through the slave-boson 
mean- field analysis of the t-J model. [ljj While the slave- 
boson theory treats strong correlation effects on average, 
it still remains a fundamental question whether the dFSD 
will be relevant to the strongly correlated regime and 
thus to the physics of high-T c cuprates. Since the dFSD 
changes the shape of the FS, possibilities of the dFSD 
are crucial to the understanding of low energy properties 
of cuprates. 

In this Letter, we study the dFSD by applying the 
exact diagonalization technique to the 2D t-J model at 
zero temperature and treat strong correlation effects rig- 
orously. We take a 4 x 4 cluster with a periodic boundary 
condition with two holes. Since the dFSD is associated 
with orientational symmetry breaking of a square lattice, 
we focus on a tendency to the orientational symmetry 
breaking. We impose a small xy-spatial anisotropy in the 
original t-J model and calculate the induced anisotropy 
of hole-hole correlations as a response. It is found that 
the response shows a sharp peak as a function of band 
parameters, which means that the system has a strong 
tendency to orientational symmetry breaking for particu- 
lar band parameters. The analysis of the momentum dis- 
tribution function indicates that this correlation develops 
when the FS is close to (tt, 0) and (0, tt). The concept of 
the dFSD, originally proposed in the slave-boson mean- 
field analysis of the t-J model^j] and a renormalization 
group analysis of the Hubbard model, [l^ is thus relevant 
even to the strongly correlated regime and in this sense 
is a generic feature of square lattice interacting electron 
systems. Possibilities of the dFSD are discussed for high- 
Tc cuprates. 

We take the 2D t-J model 

H = - t^cl a c i+Ta + ^ Jt (Si ■ S i+ r - ^n»n»- 

i, T, <7 (i,i+r) 

defined in the Fock space with no doubly occupied sites. 
The Ci„ (Si) is an electron (a spin) operator, and hi = 
So- cJcrCjcr is the number operator of electrons at site i. 
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The tr is the Zth (Z < 3) neighbor hopping integral to the 

direction t = Tj — r i; and we denote 4- — t T , tr = 

( si 

and tr — t"; the latter two are assumed isotropic. The 
J r (> 0) is the superexchange coupling between nearest- 
neighbor sites. We introduce an external anisotropy to 
t T and J T , 

t x = (l + a)t,t y = (l-a)t, (2) 
J X = {1 + -Y)J, J y = {l- 1 )J, (3) 

where we assume 7 = j^- so that the anisotropy between 
J x and J y is twice as large as that between t x and t y . As 
a response to this anisotropy, we consider an anisotropy 
of hole-hole correlation 



AC(r) = C(r) - C(r) , 
where r = (x, y) and f = (y, a;), and 



(4) 



C(r) = ^^((l-fi i )(l-ni + r)> (5) 



is a hole-hole correlation function; iV is the total number 
of lattice sites; 5(> 0) is hole density; the factor jj^gjj 
is taken to normalize J2 r C(r) to be unity. Note that 
AC(r) vanishes as long as the system has orientational 
symmetry of the lattice. 

To treat strong correlation effects rigorously, we take 
a 4 x 4 cluster with a periodic boundary condition and 
diagonalize the Hamiltonian using the Lanczos algorithm 
at zero temperature; the coordinates of the cluster are 
defined in Fig. ^a) . Although the cluster size is rather 
small for usual exact diagonalization studies, the 4x4 
cluster has square lattice symmetry and enables us to 
treat k = (tt, 0) and (0, ir). These aspects are crucial in 
the present study. We introduce two holes (S = 1/8) and 
assume that the total momentum of the wave function is 
zero. For a = 0, the system has orientational symmetry 
and hence AC(r) = 0. When a small a is introduced, 
AC(r) can become finite. It has four inequivalent points: 
(i) r = (2, 0), (ii) r = (1, 0), (3, 0), (iii) r = (2, 1), (2, 3), 
and (iv) r = (1,1), (2,2), (3,3), (3,1). HI In the latter 
case AC(r) = 0. We refer to the former three cases as 
AC m with m = 1,2, and 3, respectively. Setting a — 
0.01, we calculate AC m as a function of t" for t/J = 4 
and t'/t = -0.17. 

In Fig.^b), we see that AC m is nearly zero for \t"\ > 
0.15i, that is the hole- hole correlation is almost isotropic 
even in the presence of the external anisotropy (a/fl). 
For \t"\ < O.lt, AC m becomes noticeable and is strongly 
enhanced near t" /t — 0.05 to form a sharp peak. This 
peak is seen in all AC m , indicating that the system has a 
strong tendency to the orientational symmetry breaking 
for particular band parameters. 

In the inset of Fig. ^b), we show the following quantity 
as a function of t" to see a degree of the orientational sym- 
metry breaking: C x = XLC(r), where the summation 
is over r = (1, 0), (2, 0), (3, 0), similarly C y = £ y C(r) 
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FIG. 1: (a) The coordinates of the 4 x 4 cluster, (b) t" 
dependence of the anisotropy of hole-hole correlation AC m 
for a = 0.01, t/J = 4, and t'/t = -0.17. The inset shows 
hole-hole correlations along the x(y) direction C x ( y ) and a 
diagonal direction C xy . 



over r = (0, 1), (0, 2), (0, 3), and C xy = *£ xy C(r) over 
r = (1,1), (2, 2), (3, 3). The C x (C y ) and C xy describe 
hole-hole correlation along the horizontal (vertical) di- 
rection and a diagonal direction, respectively. Around 
t"/t = 0.05 the value of C x is about twice as large as 
that of C y . The breaking of orientational symmetry at 
t" /t 0.05 is thus very large. The hole-hole correlation, 
however, does not become one-dimensional like. The sum 
rule, ^ r ( -^( r ) = 1' suggests that the value of C xy w 0.2 
is appreciable since the mean value of C(r) is 



of C xy becomes A. 



1 and that 

Furthermore we found that C xy was 
almost the same as that for a = 0, that is the correlation 
along the diagonal direction still remained. Moving away 
from t"/t w 0.05 the difference between C x and C y gets 
smaller as expected. 

To understand the origin of the sharp peaks in 
Fig. mb), we calculate the momentum distribution func- 
tion n k = |(E CT 4 CT c feCT ) for a = [Fig. Ha)]. The 
maximal value of rik is about 0.56. This comes from 
strong correlation effects, namely the constraint that 
double occupation is forbidden at every site, which leads 
to rik < 4(1 + S). The distribution of rik is shown in 
FigS-EIb)-(d) for several choices of t"; the FS is inferred 
under the assumption of the Luttinger theorem|23j and 
may represent a typical FS in —0.2 < t" /t < —0.05 
(b), -0.05 < t"/t < 0.05 (c), and 0.05 < t"/t < 0.2 
(d), respectively. This change of FS as a funtion of 
t" agrees with one expected for noninteracting electrons 
with t, t', and t" as emphasized before . pil l25l l2*q Around 
t"/t w —0.15 the FS is expected to cross the point 
k = (ir/2, it/2). But an appreciable change of rik is 
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FIG. 2: (a) t" dependence of n k for a = 0, t/J = 4, and t' /t = 
—0.17. (b)-(d) Distribution of «k at several choices of t"; the 
k points with rik > 0.35 are considered and the diameter of 
the open circle is proportional to the value, rik — 0.35. The 
gray line is an inferred FS. (e) The change of by small 
a = 0.01 at k = (0, tt) and (tt, 0). 



not seen in Fig. HJa). This is probably due to a finite 
size effect since the total number of electron has to be 
distributed over a small number of k points constrained 
by the geometry of the cluster. A jump of n k is seen at 
t"/t = 0.05 in Fig. I^a). This jump itself is probably 
a finite size effect, which we interpret as a level cross- 
ing between a state with an electron like FS [Fig. |2Ic)] 
and a state with a hole like FS [Fig. Old)] . The free en- 
ergy shows a cusp there as a function of t". Moreover, 
the jump disappears for a^0 and nk changes smoothly, 
which will be due to realization of an "intermediate" state 
in the sense that the FS will close along the k x direction 
(electron like) and open along the k y direction (hole like) . 
To confirm that an electron like FS actually changes into 
a hole like FS near t" /t = 0.05 as inferred from Figs.^a)- 
(d), we look at the change of nk at (tt, 0) and (0, tt) by 
small a, which we define as ZVik = nk ;Q =o.oi ~ "■k ; a=o- 
In Fig.^e), we see an appreciable change of Arik around 
t"/t = 0.05 and a positive (negative) sign of An-k at (0, tt) 
[(tt, 0)]. This is expected when a FS is close to (tt, 0) 
and (0, tt) at t"/t = 0.05, where the small a moves the 
(0, tt) [(tt, 0)] point inside (outside) the FS. Therefore, 
considering that AC m shows a sharp peak at t" jt = 0.05 
[Fig-Gib)], we conclude that an electronic state with a 
FS being near (tt, 0) and (0, tt) has strong correlations 
toward orientational symmetry breaking. 

Since our anisotropy (a — 0.01) is a weak perturbation 
to the original t-J model, AC m can be a linear response 
quantity. [23] This actually holds for most of t" . When t" 
approaches 0.05i, however, the linear susceptibility grows 
rapidly and the value of a = 0.01 is no longer in the 
linear response regime; AC m begins to show saturation 



as a function of a. 

So far we have investigated the anisotropy of hole-hole 
correlation as a function of t" . A similar enhancement 
of the anisotropy would be expected as a function of t' 
(with t" — 0), since an electron like FS may cross (tt, 0) 
and (0, tt), and evolve to a hole like FS with decreasing 
t'(< 0). However, rik did not indicate a hole like FS, but 
an electron like FS [Fig. Hfc)] in a wide range —0.1 > 
t' jt > —0.43. This is probably due to larger finite size 
effects at smaller t' , since nk=(o,o) decreases drastically 
at t' jt —0.44 and becomes smaller than nk at k = 
(tt/2, 0), (tt, 0), and (tt/2, tt/2), indicating realization of 
a hole pocket around k = (0, 0). This is not expected in 
infinite systems for such parameters and in fact does not 
occur in the slave-boson mean-field analysis. Regardless 
of this possible finite size effect, AC m shows a pronounced 
enhancement near f! jt = —0.43, where we expect a FS to 
be closest to (tt, 0) and (0, tt) from the same observation 
as Fig.EJe). 

The correlation of orientational symmetry breaking 
found in Fig. IHb) is identical to the dFSD correla- 
tion that was discussed in the slave-boson mean-field 
analysis of the t-J modellll L 12(1 and in various other 
models. mBElHHlUlIiA crucial point of the 
present study is that the dFSD correlation is found in the 
exact diagonalization where strong correlation effects are 
treated exactly. Hence the <iFSD is relevant even to a 
strongly correlated regime and in this sense is a generic 
feature when a FS is close to (tt, 0) and (0, tt). 

While the present calculation has been done at a fixed 
hole density, the g?FSD correlation may be prominent in 
a certain hole doping region. The slave-boson analysis of 
the t-J model showed that the dFSD correlation was no- 
ticeable inO < 5 < 0.20 for t/J = i,t'/t = -l/6,t"/t = 
and became stronger with decreasing 5.[Tj 

The strong correlation effects are treated exactly in 
this study. But the system size is finite. Does a spon- 
taneous g?FSD take place when the system size becomes 
infinite? This is an open question. Since the g?FSD com- 
petes with other instabilities such as d-wave supercon- 
ductivity, there is a possibility that a spontaneous <£FSD 
is hindered by a more dominant instability. [TH Il9| 

Even if a spontaneous dFSD does not occur, appre- 
ciable correlations of the dFSD make a system sensitive 
to a small external anisotropy, which leads to a notice- 
able dFSD — the FS is softened. [U |U This sce- 
nario was proposed for La-based cuprates from strong 
band parameter dependence of the ciFSD correlation 
[see Fig. ITTb)]. Ill] The shape of the FS was interpreted 
in terms of the dFSD under the assumption of a cou- 
pling to the low-temperature tetragonal lattice distor- 
tion that gave a small external anisotropy to the electron 
system. [ll| Magnetic excitations were investigated on the 
basis of the dFSD scenario and turned out to catch es- 
sential features of neutron scattering data. [2^] Recently 
anisotropic neutron scattering signals were observed for 
untwinned yttrium barium copper oxides. |29| and effects 
of the g?FSD will be interesting for this material also. 
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The (iFSD correlation is the same as electronic "ne- 
matic" correlations 30J from the view of symmetry - 
orientational symmetry breaking of a square lattice. 
The possibility of nematic order has been discussed for 
high-T c cuprates. 10] It is, however, currently envisaged 
as melting of spin-charge stripe order; the underlying 
physics is different from the cZFSD. In the stripe sce- 
nario, some interaction with large momentum transfer 
near q = (tt, tt) is requisite to realize spin-charge stripe 
formation. On the other hand, the <iFSD comes from for- 
ward scattering of electrons close to the FS near (tt, 0) 
and (0, 7r); stripes are not necessary to get the nematic 
order. It should be noted that the dFSD scenario docs 
not preclude realization of stripes, which could be driven 
by interactions with large momentum transfers in the 
dFSD state. In the t-J model, however, this is unlikely. 
The hole-hole correlation is not one-dimensional even in 
the presence of the anisotropy [inset of Fig. db)]. In 
addition, various numerical studies of the t-J model re- 
ported that stripe formation was not favored [3ll l33 | es- 
pecially in the presence of t' and £",[32, 13 that is for 
realistic band parameters of cuprates. 

In summary, we have investigated a tendency to ori- 



entational symmetry breaking of the square lattice t-J 
model. To treat strong correlation effects rigorously, we 
have applied the exact diagonalization technique to a 4x4 
cluster with a periodic boundary condition with two holes 
at zero temperature. Introducing a small anisotropy to 
t and J, we have calculated the induced anisotropy of 
hole-hole correlations as a response. We have found that 
the response is strongly enhanced for particular band pa- 
rameters, meaning that the system has a strong tendency 
to orientational symmetry breaking of the square lattice. 
The analysis of momentum distribution function indi- 
cates that this correlation develops when the FS is close 
to (tt, 0) and (0, tt). The present study shows that the 
idea of a <fFSD, originally proposed in the slave-boson 
mean-field analysis of the t-J model[Tll| and a renormal- 
ization group analysis of the Hubbard model, 4.2] is rel- 
evant even to a strongly correlated regime and is thus a 
generic feature of square lattice interacting electron sys- 
tems. Possibilities of the dFSD are interesting for high-T c 
cuprates as well as other materials. |35| 
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